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Abstract Studies of the Spectral Energy Distribution of Young Stellar Objects suggest that the outer disk of FU 
Orionis objects might be self-gravitating. In this paper we propose a method to test directly whether, in these 
objects, significant deviations from Keplerian rotation occur. In a first approach, we have used a simplified model 
of the disk vertical structure that allows us to quickly bring out effects related to the disk self-gravity. We find 
that the often studied optical and near-infrared line profiles are produced too close to the central object to provide 
significant evidence for non-Keplerian rotation. Based on parameters relevant for the case of FU Ori, we show that 
high-resolution long-wavelength spectroscopy, of the far-infrared H2 pure rotational lines (sometimes observed in 
"passive" protostellar disks) and sub-mm CO lines, should be well suited to probe the rotation curve in the outer 
disk, thus measuring to what extent it is affected by the disk self-gravity. The results of the present exploratory 
paper should be extended soon to a more realistic treatment of the disk vertical structure. 

Key words, accretion, accretion disks - gravitation - stars: pre-main sequence 



1. Introduction 

Very young (Class or Class I) protostellar sources are 
thought t o be characterized by a fairly high mass accre- 
tion rate llCalvet et al.Ll200nt) . but they appear to be still 
deeply embedded in their protostellar envelopes, so that 
their disk remains hidden to direct observations. On the 
other hand, in the more evolved T Tauri stars the mass ac- 
cretion rate is generally considered to be modest, so that 
their disks are heated by the combined contribution of 
viscous dissipation and irradiation from the central star, 
making it more difficult to extract detailed information 
about t he accretion process. In this context, FU Orionis 
objects ijHartmann fc Kenvonl Il996() are a rather small 
but remarkable class of pre-main sequence stars, because 
they are ideal "laboratories" to test the process of disk ac- 
cretion during the early stages of star formation. In fact, 
the disk of FU Orionis objects, differently from that of 
Class I objects, can be studied directly from its optical 
emission, and, differently from that of T Tauri stars, is 
likely to be the site of "active" accretion. 

The distinctive feature of FU Orionis objects is that 
they undergo a violent outburst phase. During the out- 
burst, they can increase their luminosity by more than 4 
magnitudes in a period of a few years, and then slowly 
return to a quiescent state on a much longer timescale. It 
is commonly believed that such outbursts are the result of 
a significant increase of the mass accretion rate in the disk 



which usually surrounds Young Stellar Objects. Many pos- 
sible mechanisms to trigger the outburst phase have been 
discussed in the literature, am ong which are a tidal inter - 
action with a companion star jBonnell &: BastienL Il99^ . 
the onset of a thermal instability in a partially ionized 
disk, when the o uter disk is alrea dy accreting at a suffi- 
ciently high rate l|Bell et al.lll995|) . or the onset of thermal 
insta bility induced by the presence of a satellite within the 
disk llClarke fc Sveilll99fil) . The high accretion rate makes 
the emission of these systems dominated by the accretion 
luminosity, with only a minor contribution from the cen- 
tral star. 

The modeling of the Spectral Energy Distribution 
(SED) of FU Orionis objects provides an estimate of the 
product Mi,M (as discussed below in Sect. 12.111 which 
turns out to be of the order of 10~^Afg/yr. For a stellar 
mass Mi, w M©, this would correspond to an accretion 
rate as high as « lO~'*M0/yr. "Standard" models based 
on the presence of a Keplerian accretion disk lead to a 
satisfactory fit of the available photometric data only for 
wavelengths shorter than 10/im; the luminosity at longer 
wavelengths is usually much larger than expected. A flar- 
ing disk model, in which the outer disk is illuminated by 
the inner disk, could in principle explain the far-infrared 
excess, but the req uired amount of flaring turn s out to 
be often too large feenv on fc Hartmannl ll'99l[l : there- 
fore, the long- wavelength part of the SED is generally at- 
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tributed to an infalling envelope, which is heated by the 
accretion disk luminosity. 

Recently we have s hown the viability of the picture 
l)Lodato fc Bertinll200ll hereafter LB) in which the lone 
wavelength SED of FU Orionis objects is considered to 
be the signature of the effects of the disk self-gravity. 
There are already many clues that point to the impor- 
tance of the disk self-gravity in this context. Among 
these, we may recall that: (i ) Submillimeter observations 
dSandell fc Weintraubl Effll show that the disk masses in 
these systems are much higher than those of T Tauri stars 
(indicating that they are probably younger than their low- 
luminosity counterparts); (n) The modeling of th e out- 
burst in terms of a thermal instability event (Bell fc LinL 
requires a very low value for the viscosity param- 
eter a, which in turn leads to a high surface density; in- 
deed, detailed vertical structure calculations (Bell et al., 
Il997j) show that at low values of the viscosity the disk is 
marginally stable against axisymmetric gravitational dis- 
turbances already at a distance significantly smaller than 
1 AU from the central star. In fact, it has been suggested 
that FU Orionis outbursts could be triggered by a mech- 
anism that is based o n the role of the disk self-gravity 
l)Armitage et al.ll20oH) . 

In the model proposed by the far infrared excess 
is produced by two separate contributions: (i) The higher 
rotation curve of the disk, in the case where the disk is 
sufficiently massive, enhances the viscous dissipation rate; 
(ii) The extra heating required to keep the outer disk 
marginally stable with respect to gravitational instabili- 
ties makes the surface temperature higher at large radii. 

The second contribution is probably related also to 
non-local energy transport processes that should be- 
come important when global gravitational instabilities 
are present. The importance of global effects in self- 
gravitating disks ha s been brought ou t by means of nu- 
merical simulations ijRice et al.L l200,s[l . but further work 
is desired to clarify the issue of non-local transport in 
massive disks (Lodato & Rice, in preparation). Thus, for 
such second contributi on, much of the discussion f ocuses 
on theoretical aspects l|Balbus fc PapaloizouLll999|) . 

The importance of the first contribution can be 
checked by means of direct observations. In this paper 
we address the issue of detecting non-Keplerian rotation 
in protostcUar disks. In passing, we note that the non- 
Keplerian rotation observed in some AGN accretion disks 
can be successf tillv explained by the model on which this 
paper is based (|Lodato fc Bertin. .200.3) . Here we propose 
a possible measurement of the rotation of the outer disk in 
FU Orionis objects based on the analysis of mid-infrared 
and sub-millimetric spectroscopy. While sub-millimetric 
emission from FU Orio nis objects is generally c onsidered 
to come from the disk (* Weintraub et all llQQj) . alterna- 
tive scenarios are available for the interpretation of the 
far- infr ared emission. In fact , as w e have already men- 
tioned, iKenvon fc HartmannI lll 99lh attribute th e far in- 
frared emission to the envelope. Here, following ILBL we 



will assume that the far infrared emission comes from the 
disk. 

A useful diagnostics to probe the kinematical prop- 
erties of FU Orionis disks is provided by the shape of 
the observed line profiles. One merit of this kind of di- 
agnostics resides in the fact that the surface tempera- 
ture of the disk decreases with radius, and that the wave- 
length characterizing the emission from a certain annulus 
of the disk changes accordingly. Therefore, the study of 
the profiles of different lines allows us to probe the disk 
kinematics at different radial distances from the central 
accreting protostar. In particular, the observed optical 
and near-infrared liii e profiles are usually double-peaked 
I Kenvon et al I hflSSL hereafter KHH), as expected from 
a rotating disk; in addition, the peak separation decreases 
with increasing wavelength, as exp ected for a disk with 
a decreasing temperature profile ( Popham et all Il996l 
hereafter PKHN). We can thus anticipate that from the 
peak separation at very long wavelengths in the far in- 
frared and in the sub-mm, probing the outer disk, we 
should be able to tell whether the rotation is significantly 
different from Keplerian, hence testing the predictions of 
a model in which the disk self-gravity affects the shape of 
the rotation curve. 

This kind of observations will require high spectral 
resolution (in order to detect small velocity differences) 
and high spatial resolution (in order to disentangle dif- 
ferent kinematical comnponents). In fact, many protostel- 
lar disks are known to possess strong molecular outflows. 
Therefore, a possible source of confusion in the shape of 
the line profiles may come from the kinematics of the out- 
flows, that could result in asymmetries in the shape of 
the line profile. In this respect, the high spatial resolu- 
tion that can be achieyed w ith interferometric techniques 
dOnillotean fc nutrevLlT99^ . and, in the near future, with 
ALMA will eventually lead to clear-cut tests on the role of 
disk self-gravity in the rotation curve of protostellar disks. 

The paper is organized as follows: in Section |21 we de- 
scribe our method to test the rotation properties of the 
outer accretion disk; in Section |3| we describe the ba- 
sic properties of the adopted accretion disk model based 
on parameters relevant for FU Ori; in Section 0] we de- 
scribe our expectations for optical and NIR line profiles; 
in Section [S] we consider the case of mid-infrared pure 
rotational H2 line profiles; in Section we describe the 
expected shape of sub-mm CO line emission; in Section d 
we draw our conclusions. 

2. Measuring deviations from Keplerian rotation 
in spatially unresolved protostellar disks 

The method presented in this Section is an extension of 
the method described bv IKHhI . Our main assumption is 
that the SED of FU Orionis objects is produced by an 
optically thick accretion disk up to wavelengths of the 
order of 100/im, consistent w"ith studies of the millim etric 
continuum in these systems ijWeintraub et allll99ll) . The 
basic procedure can then be summarized as follows: 
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1. The surface temperature profile Ts(r) is obtained from 
a parametric fit to the SED; note that in this paper 
we are not interested in explaining the physical origin 
of this profile, which we addressed in a previous article 
(LB). Given a value of the inclination angle i, this fit 
leads to a determination of the inner and outer radii of 
the disk, ri„ and Touti and of the product MM^,, where 
M is the mass accretion rate and is the mass of 
the central object. 

2. The profiles of the optical-NIR lines (such as some CO 
absorption bands, seen in FU Orionis objects), which 
are produced mostly in the inner disk, are used to de- 
termine the rotational velocity in the inner Keplerian 
disk. For a given value of the inclination z, fitting the 
observed line shapes yields a measurement of Af^/ri„. 
As a result of the first two steps, we obtain separate 
estimates of and Af . So far, the procedure is very 
close to that of lKHHl 

3. In significantly massive disks, the line profiles at longer 
wavelengths should be broader than expected from 
a Keplerian extrapolation of the rotation curve mea- 
sured in the previous step. Fitting the observed mid- 
infrared and sub-millimeter line pro files with the self- 

Jravitating disk model described in iBertin fc LodatJ 
1999. hereafter BL) will then provide a value for the 
lengthscale that marks the transition to the non- 
Keplerian rotation regime (see Subsection 12 . 21 belowl . 
This fit will test directly the importance of the disk 
self-gravity, independently of the arguments that can 
be made in order to explain the infrared excess in the 
SED. Furthermore, since depends on M*, M, and 
on the viscosity parameter a (see Eq. © below), this 
fit will also lead to a measurement of the viscosity pa- 
rameter a. 

2.1. Determining the radial temperature profile 

We start by introducing the dimensionless radial coordi- 
nate X = r/rin and the temperature scale Tq from the 
Keplerian model temperature profile: 



^bTKt) = 



3 GM^M 
8^ rl 



(1) 



where ctb is the Stefan-Boltzmann constant. This simple 
temperature profile is able to reproduce the correct spec- 
tral index of the SED at wavelengths smaller than lO/im 
and thus is adequate at small radii. 

In the following, we would like to test the self- 
gravitating disk hypothesis by studying directly the rota- 
tion properties of the disk, independently of a discussion 
of the detailed physical processes that justify the observed 
SED. Therefore we consider a parametric description of 
the temperature profile, in l ine with many other invest iga- 
tions of YSO disks (e.g., see iBeckwith fc Sargenll l993' and 
lOsterloh fc Beckwithlll99fl who studied the sub-mm con- 
tinuum of YSOs, or iThi et aljEnoH who studied emission 
lines). For the purpose, we thus consider a simple com- 
bination of two power laws, one describing the inner disk 



(and with temperature therefore proportional to r ^/■*), 
and the other describing the outer disk, allowing for devi- 



ations from the r 
Tt{r)^ne\x), 
where 



"^/^ dependence: 



Hx) = i 



1+ - 



(2) 



(3) 



Xt = Tt/rin being a dimensionless radius, beyond which 
the temperature profile departs from the "standard" pro- 
file, and (3 an appropriate power-law index. 

Under the assumption that the disk emission is opti- 
cally thick, the SED can be constructed by integrating 
in radius the blackbody spectrum. By analogy with what 
described in lLBl we then introduce a reference frequency 

— kTo/h and rescale the frequency to v = v/vq. For a 
given disk inclination i, the SED is thus given by: 



xda 



where 



in 



120 2 4 

— g-r^crsTo cos I, 



(4) 



(5) 



Do is the distance to the YSO, and Xout — ^ out /Tin- 

A fit to the SED of observed objects will therefore de- 
termine the fit parameters Xouti To, Lq, /?, and Xt. The 
combined knowledge of Lq and Tg (assuming the inclina- 
tion angle of the disk to be known) allows us to obtain Tin 
and the product M^M (see Eqs. (QJ and (0)). 

2.2. The rotation curve in the presence of a 
self-gravitating disk 

We now consider the contribution of the di sk se lf-gravity 
to the rotation curve. Based on the results of lBll we intro- 
duce a typical lengthscale , the radius beyond which de- 
viations from Keplerian rotation become significant, and 
a typical velocity scale Vo, defined as follows: 



CM 



-2/3 



GM 



2/3 



(6) 



(7) 



Here Q is a p arameter that, for simplicity, will be set to 
unity (s eelBlJl. The quantity a is t he parameter that de- 
fines the lShakura fc SunvaevI l)l973|) viscosity prescription. 
The rotation curve of self-regulated, self-gravitating accre- 
tion disks is given by: 



GAL 



- + —{l + $^ix/Xs)) 
X Xg 



where = G Af^./r^n, Xg = rg/rm and 0^ is defined in 
Eq. (4) oflml 



(8) 
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In the following, we will sometimes use the short, sim- 
plified expression "self-gravitating rotation curve" (as op- 
posed to Keplerian rotation curve) when we will refer to 
the rotation curve derived from Eq. ||SJ) for a model that 
includes the effects of the disk self-gravity. 

2.3. The shape of the line profiles 

In the present analysis we focus on the line shape. 
Therefore, we will consider only the relative contribution 
of the different annuli of the disk to the line intensity. 

We will express the line profiles as a function of the 
velocity shift Aw, corresponding to the frequency shift 
Av from the line center vi, so that Av/c = Av/ui. The 
Doppler line profile (f>{Av, r) at a given radius of a disk 
rotating with velocity V{r) is given by: 



4>{Av, r) 



ff (Av)^ < V^lir) 



otherwise 



where Vios{r) = V{r) sini is the maximum line-of-sight ro- 
tational velocity. Obviously, this expression takes into ac- 
count only the effects related to the disk ordered motions; 
intrinsic effects associated with the microscopic physics of 
line emission and with the thermal motions inside the disk 
are not included. 

We convolve this line profile with a Gaussian profile 
Q-{^v/vb) ^ either to simulate the effects of a finite in- 
strumental spectral resolution i?, so that Vb/c — 1/R, or 
to describe turbulent or thermal motions of the emitting 
species. The global line profile will be given by integrating 
the convolved line profile (f>{Av,r) over radius, weighted 
by an appropriate weight function l{r, vi) which expresses 
the relative contribution of the different annuli of the disk 
to the line emission: 



F{Av) 



l{r, i'i)(j){Av, r)27rrdr 



(9) 



The function Z(r, vi) is specified once the details of the 
emission or of the absorption are known. At wavelengths 
where the continuum emission is optically thin, we expect 
to find emission lines; on the other hand, where the con- 
tinuum is optically thick, emission lines are produced only 
if a temperature inversion in the upper layers of the disk is 
present, due, for e xample, to the effect of irradiation from 
the central object l)Calvet et al.lflo'Qlbl) . If we consider FU 
Orionis objects, since the accretion rate is high and irra- 
diation is not expected to be efficient, we thus expect to 
have absorption lines at optical and infrared wavelengths 
l)Ca,1vet et a,1.Lll ah . Near infrared CO absorption bands 
from th e disk are indeed a typical feature in FU Orionis 
objects However, irradiation from the inner disk 

may give rise to a hot atmosphere above the disk, where 
emission lines can be produced. Therefore, in some cases 
(for example, for the optically thin FIR H2 molecular lines, 
see SectionEJ we have considered both cases of absorption 
and emission. 



In this simplified analysis we will consider for emission 
lines the two limiting cases of optically thick (which can be 
appropriate for sub-millimetric CO rotational lines) and 
optically thin emission (best suited for FIR pure rotational 
H2 emission lines, see below), so as to keep the radiative 
transport as simple as possible. 

1. Optically thick emission lines. In this case 
the brightness temperature of the emission line co- 
incides with the physic al temperature Ts{r) (see 
iBeckwith fc Sargentlll99.'lh so that we can assume that 
the different annuli of the disk emit proportionally to 
their contribution to blackbody emission at the wave- 
length of the line center: 



l{r,vi) cx 



1 



1 



(10) 



The expected line profile can thus be obtained directly, 
based on the results of modeling of the SED, as de- 
scribed in Sec. 12.11 
2. Optically thin emission lines. In this case the rel- 
ative contribution of the different annuli requires one 
additional piece of information, about the surface den- 
sity radial profile <j(r). Apart from numerical factors 
that do not depend on the radial coordinate or on the 
temperature, we can express /(r, z/j) as: 



/(r, vi) cx a{r)- 



~E^p/kT,{r) 

Q{Ts{r)) ■ 



(11) 



where E^p is the energy of the upper level of the transi- 
tion, and Q{Ts{r)) is the partition function of the line 
emitting specie s (for example the H2 molecule; see also 
iThi et alJ200ll) . For simplicity, in this formula we have 
made the assumption that the line emitting species is 
at the same temperature as the matter responsible for 
the continuum emission (i.e. mostly the dust in the 
upper layers of the disk). However, it should be noted 
that this assumption would lead to very weak emis- 
sion lines, difficult to observe. This has been partially 
confirmed by FIR observations of pure rotational il ? 
emission lines in quiescent YSOs l|Richter et al.ll2002j) . 
that indeed reveal that these emission lines, if present 
at all, have a very weak line strength. 
A more detailed discussion of the radiative transfer 
would require a complete accretion disk model. In the 
present simplified investigation, in order to proceed 
further from Eq. Hll(l . we may consider a surface den- 
sity profile of the form a r~'^ . The power law index 
7 can be derived from spatially resolved observations 
of disks at mm wavelengths, where the disk emission 
is optically thin (Guillo teau fc Dutre v*. 'l998'). and it 
turns out to be of ten close to unity ifTesti et al.. 2003 
IPietu et all l2003l) (even for evo lved objects, such as 
TW Hya; see lWilner et al.ll2000|) . Therefore, the case 
7=1 will be our reference case (see also the discussion 
is Section lsT^ below) . [Interestingly, the self-gravitating 
disk solution of BL indeed predicts 7 = 1, but, as we 
did when we introduced the properties of the SED, in 
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this paper we do not wish to address the problem of 
interpreting the observed profile.] 
3. Absorption lines. In this case we follow IPKHnI and 
make the simplifying assumption that the equivalent 
width of the line does not change with Tg, so that 
/(r, i/j) is proportional to the blackbody emission and 
will be therefore describe d by Eg . H10|l . as for optically 
thick emission lines. As IPKHNI note, this approach 
may not be valid above some threshold temperature 
(dependent on the absorption line), where the relevant 
molecules producing absorption disappear. We have 
therefore set the equivalent width to zero above this 
temperature (see also Sec. 0] and |5l below). 
Other models of the absorption lines in FU Orionis 
include the use of stellar template photospheres with 
different temperatures, simulating the disk radial tem- 
perature profile (the disk is divided into rings, and 
with each ring a stellar template spectrum is associ- 
ated, characterized by effective temperature e qual t o 
the disk surface temperature at that radius; IKHHI ). 
or of the detailed sol ution of the radiative transfer in 
the disk atmosphere l) Gal vet et al.l Il991al) . However, 
such refined calculations show only minor differences 
with respect to estimates based on the simplest as- 
sumptions, at least for optical and near-infrared lines. 
More caution should be taken when applying the sim- 
ple constant equivalent width assumption to longer 
wavelength lines (such as the H2 lines that we are go- 
ing to consider in SeclSJ. In this first analysis we have 
decided to follow the simplified approach, being aware 
that for the final comparison with observations more 
realistic assumptions about the radiative transfer in 
the disk should be considered. 



2.4. Other factors expected to affect the shape of the 
line profiles 

Before concluding this Section, a cautionary remark is 
in order. The main objective of the present paper is to 
show quantitatively what kind of spectroscopic observa- 
tions could carry significant evidence for the presence of 
non-Keplerian motions in protostellar disks. In view of this 
goal, we have decided to consider a simplified description 
of the vertical structure and of the radiative transfer in the 
disk. Obviously, for the subsequent goal of interpreting the 
observations of some specific objects, many other effects 
concerning the dynamics, the physics and the chemistry 
of protostellar disks should be considered carefully. These 
include the possible interplay between rotat ion and infall 
in d etermini ng the s hape of the line profile ijllogerheiideL 
I2OOI: Belloche et al.[ 12002'!. th e effect of molecular deple- 
tion fsee lMarkwick et alJl20oi and comment at the end of 
Section O, and the effect of other sources of gravitational 
field in disturbing the rotation curve. This final issue can 
be particularly important in the case of FU Orionis ob- 
jects, for which the tidal interaction wi th a close "compan- 
ion" has been sometimes suggested (^Bonnell fc BastienL 



as a possi ble way to trigger the outbu rst. However, 
event statistics llHartmann fc KenvonLll99^ suggest that 
FU Orionis might be a repetitive phenomenon, with a re- 
currence period of « 10'* — 10^ yrs. The "companion" 
should therefore lie in a very eccentric orbit, so that it 
should cause only a minor disturbance to the shape of the 
line profiles after the outburst has been triggered. 

3. Basic parameters for a model suggested by the 
SED of FU Ori 

In this Section we apply the initial step of the proce- 
dure outlined above to th e SE P of FU Ori, the proto- 
typical FU Orionis object. ILBI have already modeled the 
SED of this object in terms of a self-gravitating accretion 
disk, obtaining a satisfactory fit to the available data up 
to lOO^m. Therefore, this object is a good candidate for 
testing its rotation characteristics and possible deviations 
from Keplerian rotation. In Table H we report the basic 
physical parameters for FU Ori obtained from the mod- 
efing of its SED in previous analyses (i.e. IKHHL IPKHNL 
and and from the parametric modeling described in 
Subsect. 2.1 and thus adopted in this work (see below). 

The different models cannot be easily compared to 
each other, because of the different parameters assumed 
(distance, extinction and inclination angle) and of the dif- 
ferent modeling considered for the emission in the disk 
boundary layer. In particular, the main characteristics of 
the models are: 

- KHH assume that the net angular momentum fiux 
in the disk is not vanishing and is given by the re- 
quirement of zero torque at the inner disk radius. The 
boundary layer is assumed to be at a constant temper- 
ature, corresponding to the maximum disk tempera- 
ture. 

- PKHN assume that the net angular momentum fiux in 
the disk is vanishing and include the boundary layer 
emission, treating it as a low-radiative-efficiency fiow. 

- LB (and this work) assume that the net angular mo- 
mentum flux in the disk is vanishing, but do not in- 
clude the boundary layer emission. 

The differences in the choice of the visual extinction 
Ay (in magnitudes) and of the distance Dq only induce 
slight changes in the estimate of the product M^^M , when 
similar values for the inclination are assumed. On the 
other hand, the different treatment of the inner disk emis- 
sion leads to a different estimate of the inner disk radius 
in the three models, so that r™^ < rf^^'^ < . We 
interpret the relatively large value of the inner radius es- 
timated by LB and in this work as a consequence of the 
fact that, with respect to PKHN, we did not include the 
contribution of the radiatively inefficient boundary layer, 
so that our inner disk radius is to be considered closer 
to the transition from the thin disk solution to the hot 
boundary layer, rather than to the actual stellar radius. 
Another important difference between KHH and PKHN 
concerns the central mass estimate: the boundary layer in 
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Table 1. Accretion disk parameters for FU Ori, for the disk models of KHH, PKHN, LB, and this work. Values 
marked with an asterisk are assumed. 







M/(M0/yr) 
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0.34 
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PKHN 


0.7 
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LB 
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0.82 10"* 


0.82 10"* 
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0.65 


2* 


550* 


this work 


0.7* 


1.4 10"* 


1 10"* 


11 


0.5* 


2.2* 


500* 
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Log((//hz) 

Figure 1. Spectral Energy Distribution of FU Ori. Th e 
triangles show the data (from lKenvon fc Hartmannll99ll) , 
while the curves show the best fit model obtained in this 
work (solid line) and by LB (dotted line). 



the PKHN model rotates slower than Keplerian, due to the 
effect of the increased pressure gradient in the boundary 
layer. Thus, to obtain the same rotational velocity in the 
inner disk (to match the optical line profiles; see below) a 
larger central mass is required. 

We will follow the work of PKHN as much as possible. 
Therefore, we start by assuming the extinction coefficient, 
distance and inclination angle provided in their analysis. 
In particular, the inclination angle assumed below is such 
that cos i = 0.5, slightly smaller than the value cos i = 0.65 
derived by LB. 

With the parametric model described in Section 12.11 
we have fitted the avail able photometric data (from 
iKenvon fc HartmanrJl99ll), corrected by m eans of a stan- 
dard dereddening function l|Cardelli et al l lT989ft . The de- 
rived values for the product M^,M and are shown in 
Table U] Figure [T] shows the observed SED of FU Ori to- 
gether with our best fit (solid line) and the best fit ob- 
tained by LB with a self-gravitating disk model (dotted 
line). The best fit values for xt and f3 (see Section l^T^ are: 
Xt ~ 261 (which corresponds to rt « 13 AU) and /? = 3. 



^4 



r/AU 

Figure 2. Surface temperature profile of FU Ori obtained 
in this work from a parametric model of the SED (solid 
line), compared to that obtained by LB (dotted line). 

In view of the goals of this paper, of the inhomogeneity 
of the SED data, and especially of the detailed discussion 
provided in KHH, PKHN, and LB, we will not comment 
here on the issue of setting uncertainties in the value of 
the derived parameters recorded in Table 1 . We only note 
that the relative poor quality of the SED data might re- 
sult is some ambiguities in the values of the estima ted 
parameters, as already noted bv lThamm et al.l l(l994l) . 

Note that (3 — 3 implies that the surface temperature 
profile is flat in the outer disk. In Fig.|21the temperature 
profile of our best fit parametric model (solid line) is dis- 
played together with the best fit model profile obtained 
by LB (dotted fine). 

4. Optical-NIR line profiles 

The contribution of the disk mass to the rotation curve de- 
pends on the size of the transition radius (see Eq. ©) 
in relation to the region of the disk actually probed by the 
available spectroscopy. In this Section we will show that, 
in the case of a model suggested by the study of the SED 
of FU Ori, rg is large enough so that the disk self-gravity 
has little or no influence on the optical-NIR line profiles. 
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Figures. Optical (6170A) and near-IR (2.2/im) predicted 
line profiles for FU Orionis, in the Keplerian (solid line) 
and self-gravitating (dotted line) model, with rg « 18 AU. 
The curves practically overlap and no differences between 
the two descriptions are noted. 

which are mostly produced in the inner disk. In practice, 
from the peak separation of the observed line profiles at 
these relatively short wavelengths, one can measure the 
ratio Mi, /rim thus preparing the ground for a firm ex- 
pectation of effects related to deviations from Keplerian 
rotation in spectroscopic studies at larger wavelengths (see 
following Section). This Section thus addresses step 2 of 
the procedure outlined in Sec. |2] 

Here we follow PKHN and assume their value for the 
central mass, Af* = O.TMq (PKHN found significantly 
larger values of Af^, with respect to KHH who did not 
include a detailed model of the inner boundary layer). 
Assuming a — 0.034 (the best fit value obtained by LB), 
we would obtain « 350ri„ w 18 AU. Actually, as a 
result of the different estimates of M^, and M considered 
here with respect to LB, this value for is smaller than 
that obtained by LB, leading to a higher and possibly 
unrealistic total disk mass (w 2Af+). We have thus also 
considered the case in which a — 0.1, leading to « 
700ri„ « 36 AU, better in line with the best model of LB. 

We refer to the same (absorption) lines as in PKHN 
and KHH, i.e. metal lines at 6170 A(at a resolution 
R = 24000) and CO lines at 2.2/im (at a resolution 
R = 20000; the spectral resolutions are taken to be the 
same as in PKHN). The threshold temperatures above 
which the equivalent width of the line is taken to be van- 
ishing are 8000K for the optical lines, and 5000K for the 
near-infrared lines. In Fig.|3]we plot, for w 18 AU, both 
the profile resulting from strictly Keplerian rotation (solid 
line) and the one resulting from the self-gravitating model 



(dotted line). The two different profiles are clearly indis- 
tinguishable. Note that the quantity usually referred to is 
not the actual absorption line profile, but the peak of the 
cross-correlation function of many absorption lines in the 
same wavelength range. This quantity is preferred from an 
observational point of view because it gives a way to ob- 
tain higher signal-to-noise ratios; physically, it represents 
the average shape of the line profile in a small wavelength 
range. 

A similar plot for the case in which w 36 AU would 
again show, as expected, no significant differences in the 
shape of the line profiles between the strictly Keplerian 
model and the self-gravitating disk model. 

5. Probing the outer disk rotation with H2 line 
profiles 

In order to probe properties of the outer disk, so as to 
have a chance for detecting significant deviations from 
Keplerian rotation, one should consider line profiles at 
longer wavelengths (see Fig.|21l. 

Molecular hydrogen pure rotational emission lines have 
been observed with the Infrare d Space Observatory (ISO 
in many protos tellar systems l|van Dishoeck et"ail Il99 
iThi et aL[l200li) . However, due to the low angular resolu- 
tion of ISO, it is not clear whether the H2 emission can 
be actually associated with the disk. Indeed, these ISO 
obser vations were not con firmed at higher spatial resolu- 
tion ^Richter et al.U2002() . The weakness of Ht emission 
lines could indicate that spatial and temperature sepa- 
rations between the gas and the dust in the disk are 
small. Molecular hydrogen is a homonuclear molecule, so 
that the rotational emission lines result from the electric 
quadrupole. Therefore, they have a very low Einstein A 
coefficient, which makes them optically thin up to high 
column densities (of the order of lO^^cm^^). We will there- 
fore assume that these lines are optically thin, although we 
should be aware that, if the self-gravitating disk hypothe- 
sis is correct, then the surface density of the disk is going 
to be high as well, so that the optically thin assumption 
should eventually be re-examined. The basic line param- 
eters of the H2 rotational lines are summarized in Table 
11 

A complete analysis would require the discussion of 
the full disk model, with an assumed vertical temperature 
profile, and the solution of the radiative transfer equation 
in the vertical direction, taking into account both viscous 
dissi pation and exte rnal irradiation as heating sources 
( Calvet et al l . Il991al) . Here, for simplicity, we will only 
consider the two extreme cases of absorption lines (appro- 
priate for the case in which viscous dissipation dominates) 
and of emission lines (that may be produced as an effect 
of irradiation of the outer disk by the inner disk), under 
the hypothesis that the line emitting gas is at the same 
temperature as the continuum emitting dust, i.e. that the 
gas temperature is given by Ts as derived from the model- 
ing of the SED. The line profile will be obtained following 
the procedure described in Section lOI 
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Table 2. Pure rot ational H? line p arameters: central wavelength, energy of the upper level (in Kelvin), Einstein 
A-coefficient. From lThi et~alT ()200l|i . 
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Figure 4. Molecular hydrogen absorption line profiles for disk models with self-gravitating (with ~ 18 AU, solid 
line) and Keplerian (dotted line) rotation curve. The upper panel shows the S(0) line at 28/im, with a resolution of 
5km/sec (left) and 4km/sec (right); the lower panel shows the S(l) line at 17/im, with a resolution of 5km/sec (left) 
and 4km/sec (right). 



5.1. Absorption lines 

Figure 0] shows, for ~ 18 AU, the predicted line profiles 
for the Keplerian and non-Keplcrian rotation curves, for 
the S(0) and S(l) absorption lines (here plotted "upside- 
down" to allow an immediate comparison with the optical 
and NIR peaks), at different spectral resolutions. We have 



used the same input parameters as for the optical and 
near-infrared line profiles. The shape of the absorption 
line profiles is not strongly affected by the choice of the 
threshold temperature above which the equivalent width 
is taken to vanish. Lowering the threshold from lOOOOK 
to lOOOK only leads to a slight effect on the wings of the 
lines. This is because the inner regions of the disk give 
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Figure 5. Same as Fig. 0] but with « 36 AU. 

only a small contribution to the line profile, due to the 
small area they cover. 



In Fig. Owe show the same line profiles, but for the 
larger value of r^, i.e. w 36 AU. From these two figures, 
it is clear that, in order to detect the small changes in the 
rotation curve associated with the effects of the disk self- 
gravity, we should reach very high spectral resolution (see 
Section 15.31 for a discussion of the feasibility of such ob- 
servations). Another interesting point is that the S(0) and 
S(l) line profiles, at the resolution of 5 and 4 km/sec, can 
be single peaked in the strictly Keplerian case and dou- 
ble peaked in the self-gravitating case. This would clearly 
distinguish between the two different pictures, if an ade- 
quate signal-to-noise ratio could be achieved. As expected, 
when a larger value of is assumed, the effects of the disk 
self-gravity are less prominent. 



T 
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5.2. Emission lines 

Figure shows the predicted emission line profiles for 
the Keplerian and self-gravitating rotation curves, for the 
S(0), S(l) and S(2) lines at different spectral resolutions 
for the conservative case in which Vs ~ 36 AU. 

The general features of the shape of the lines are sim- 
ilar to the case of absorption: the lines are broader in 
the self-gravitating case than the corresponding lines in 
Keplerian disks, and they are often double-peaked, also 
in cases where the Keplerian lines are single-peaked. In 
the model where the transition radius is smaller (not 
shown here), the effects of self-gravity are more prominent, 
as expected. 

The shape of the S(2) line profile is practically not in- 
fiuenced by the modification of the rotation curve (i.e. it 
is mostly produced in a region where Keplerian rotation 
dominates). Based on the analysis of this and of the pre- 
vious Section, we conclude that the best lines to probe de- 
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viations from Keplerian rotation would be S(0) and S(l), 
because they would display some difference between the 
two competing models, with peak separation large enough 
to be detected with reasonable spectral resolution. 

5.3. Feasibility 

The above results obtained here show how deviations 
from Keplerian rotation can be detected with the anal- 
ysis of the shape of the molecular hydrogen line pro- 
files in the mid-infrared, provided that we can achieve 
high spectral resolution and a high signal-to-noise ratio. 
The Texas Echelon Cross Echelle Spectrograph (TEXES) 
operates in the wavelength range from 5.5 to 28.5 fim 
and can reach a resolution as high as 2-3 km/sec at the 
shortest wavelengths and of w 5 km/sec at 17/im. It 
has been mounted on the NASA's 3m Infrared Telescope 
Facility, to confirm the existence of molecular hydrogen 
rotational lines apparently observed with ISO in low- 
mass, non-outbu rsting objects, such as some T Tauri and 
Herbig Ac stars ijRichter et al.Ll2002l) . In addition, EXES, 
the companion instrument of TEXES, is expected to be 
placed on SOFIA (Stratospheric Observatory for Infrared 
Astronomy). Here we have demonstrated that, in princi- 
ple, these two instruments have sufficient spectral resolu- 
tion to measure this effect in FU Orionis disks. 

If we consider instead the capabilities of ISO, the maxi- 
mum spectral resolution achievable is 10 km/sec. The pre- 
dicted S(0) line profiles in FU Ori with the lower resolu- 
tion of 10 km/sec is no longer double-peaked. In the self- 
gravitating case there is only a slight broadening of the 
line profile with respect to the Keplerian case. We thus 
consider it unlikely that the signature of non-Keplerian 
rotation can be found in ISO observations. 

Actually, the most serious problem with these mid- 
infrared line observations is to achieve a sufficiently high 
signal-to-noise ratio. For the case of emission lines, if 
the line emitting gas and the continuum emitting mat- 
ter are approximately at the same temperature, we expect 
only weak emission lines to be present (as confirmed by 
iRichter et al.ll200^ . For the case of absorption, the small 
absorption coefficient of the pure rotational H2 transitions 
also leads to weak absorption lines. 

5.4. Robustness of the H2 emission lines results 

The optical thickness of the relevant lines depends on the 
surface density profile, which cannot be directly inferred 
from the modeling of the SED alone. Observations of pro- 
tostellar disks, consistent with our self-gravitating disk 
model (BL), suggest a disk density profile of the form 
a oc r^'^ , with 7 « 1; in the previous Sections, we have 
indeed set the surface density power law index to be 7 = 1. 

The H2 rotational emission line are optically thin up 
to relatively high densities, due to their low Einstein's A- 
coefficient. FU Orionis objects are likely to be younger, 
and denser, than T Tauri and protoplanetary disks, so 



that for them we may even have to consider the possibility 
that the H2 emission is optically thick. In this Section 
we discuss the robustness of the results presented in the 
previous Section, with respect to the uncertainties on the 
optical thickness of the disk and on its surface density 
profile. 

In Fig. [7| the local emissivity of every annulus of the 
disk, for the S(0) and S(l) fines, is plotted as a function of 
radius for three different cases, i.e.: (i) optically thin emis- 
sion, 7=1, (zi) optically thin emission, 7 = 1.5, and (m) 
optically thick emission. When the steepness of the sur- 
face density profile is increased, the peak of the emission 
for a certain line occurs at relatively smaller radii, as a 
result of the increased surface density. On the other hand, 
optically thick emission tends to contribute more than op- 
tically thin emission at larger radii, because in our models 
the temperature profile in the outer disk is approximately 
constant, while the optically thin emissivity rapidly de- 
creases, being proportional to the surface density profile. 

In Fig. IHl we show how the S(l) line profile changes 
when a different slope of the disk density is assumed. A 
steeper density profile (here we consider 7 — 1.5) results 
in a larger peak separation. This is due to the fact that 
a steeper surface density profile lets the inner parts of 
the disk, which rotate faster, contribute more prominently, 
with respect to the case of a less steep profile. This identi- 
fies a potential problem in the observational test that we 
are proposing here (for the case of optically thin line pro- 
files). In fact, the increased contribution of the inner disk 
to the line emission when a steeper surface density profile 
is assumed may be such that the S(l) line is produced 
mostly at radii where the rotation curve is Keplerian. The 
resulting modifications to the line profiles due to the disk 
self-gravity would therefore become only barely visible in 
this case. 

In practice, there is a sort of degeneracy in the ap- 
pearance of optically thin line profiles with respect to the 
assumed model: a broader line profile can be traced either 
to the effects of the disk self-gravity, which enhances the 
rotation at large radii, or to the steepness of the disk den- 
sity profile, which makes it possible for a certain emission 
line to originate from smaller radii, where the rotation 
velocity is larger. 

6. Sub-mm CO line profiles 

The discussion of the previous Section shows that some 
H2 emission lines can be a useful tool to probe the rota- 
tion of the outer disk, although confusion may arise in the 
interpretation of optically thin lines because, under cer- 
tain circumstances, deviations from the typical a ^ 1/r 
disk density profile can mimick those associated with de- 
viations from Keplerian rotation. 

To derive more convincing conclusions about the rota- 
tion of the outer disk, it would thus be important to con- 
sider long- wavelength lines that are optically thick, so that 
the additional assumptions concerning the surface density 
profile are not needed. In this respect, some sub-mm CO 
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lines are very inte resting, because they are expected to 
be op tically thick l)Beckwith fc Sargentl Il993t IXhi et all 
l200ltl and can be studied with yery high spectral resolu- 
tion. These global CO line profiles, observed in T Tauri 
and protoplanetary disks by a number of authors, usually 
show the typical do uble-peaked profile ijXhi et all I2OOII: 
iDutrev et allll997(l . 

Interestingly, millimctric CO line emission has als o 
been observed in V1331 Cyg l|McMuldroch et all Il99l . 
an "FU Orionis object between outbursts", which shares 
many features with FU Orionis objects, but has a smaller 
luminosity. The estimated disk mass in V1331 Cyg is 
PS 0.5Mq, a further clue to the fact that FU Orionis disk 
mass may be fairly high. 

As an example, we consider the ^^CO 1-0 at 110 GHz, 
which is commonly observed in protostellar disks. We as- 
sume that the main limiting factor for the spectral res- 
olution is given by turbulent or thermal broadening and 
we take the factor Vb in the Gaussian with which the line 
profile is convolved (see Section to be 1 km/sec. The 
resulting line profiles, for the two different choices of r^, 
are shown in Fig. ^ The effect of self-gravity is clear in 
both cases. 

In any case, as noted in the Introduction, care should 
be taken (especially for younger objects, such as FU Ori) 
in order to properly identify the various kinematical com- 
ponents (such as winds and outflows) that might con- 
tribute to the shape of the line profile. Another source 
of uncertainty related to CO lines is that of molecular 
depletion, resulting from the freezing out of molecular 
species onto dust grains. Many studies (see, for example, 
[ Markwic k ct al. 2002) have shown that, in the outer and 
colder disk, the column density of many species can be re- 
duced significantly. However, these studies refer mostly to 
parameters appropriate for T Tauri disks, while we expect 
that FU Orionis disks, which are characterized by a much 
higher temperature, will be less subject to depletion. 

7. Conclusions 

The study of the line profiles produced in Young Stellar 
Objects has traditionally played a very important role in 
determining the physical conditions of circumstellar mat- 
ter. In particular, it has been used to confirm the accretion 
disk origin of the luminosity of FU Orionis objects (KHIJ) 
and to test the Kepler i an ro tation in many T Tauri stars 
ifGuilloteau fc Dutrevl Il994|) . In addition, spectroscopic 
observations may be able t o probe the vertica l structure of 
the disk (see, for example. ICheng et al.lll98^ for the disks 
in dwarf novae). 

Spatially resolved sub-mm interferometric observa- 
tions of T Tauri stars have already provided inter- 
esting information about the outer disk k inematics 
ijGuilloteau fc Dutreyi. .l998; Simon et alll200l|) . Such ra- 
dio observations can reach an angular resolution of « 1". 
However, for objects such as FU Ori, which lies at a dis- 
tance of ^ 500 pc, the angular size of the disk (which 



should extend out to 50-60 AU, judging from the model- 
ing of the SED) is likely to be as small as « 0.1". 

In this paper we have exploited the wavelength de- 
pendence of the radial position of the peak of the line 
emissivity of the disk, to obtain additional information on 
the kinematical properties of the outer disk in FU Orionis 
objects, based on the shape of global line profiles. Our 
main purpose has been to define a clear cut test about 
the importance of the disk self-gravity in determining the 
rotation curve of the accretion disk in the early stages of 
star formation. 

The method presented in this paper is based on the 
complementary information that can be obtained from 
the analysis of the optical and near-infrared line profiles 
(which probe the inner rotation curve, and hence set the 
value of the mass of the central object) and the analysis of 
high spectral resolution line profiles at long-wavelengths 
(which probe the rotation curve in the oute r disk). This is 
a generaliz ation of the method d iscussed bv lKenvon et alJ 
(jl988.) and lPopham et al.l l|l99(j) . 

The high spectral resolution that can be obtained 
with the ground based and already operating TEXES in- 
strument, and with its air-borne companion instrument 
EXES that wi l l be o n board of the SOFIA observatory 
l|Richter et all I2OO2I) . may indeed give us the way to 
distinguish between Keplerian and non-Keplerian rota- 
tion based on the analysis of pure rotational H2 lines. 
This in an important test about the viability of the self- 
gravitating disk model in FU Orionis obj ects, independent 
of the results of the SED modeling bv iLodato fc Bertii] 
l)200ll) . In addition, once we accept the self-gravitating 
disk model, this test will give us an alternative, indepen- 
dent way to measure the viscosity parameter a, thus pro- 
viding insights into the physical mechanisms that operate 
in protostellar disks. 

As far as the shape of optically thin line profiles is 
concerned, we have also considered the robustness of our 
conclusions with respect to changes in the surface density 
profile of the disk. We have found that changes in the 
surface density profile could in principle mimick the effects 
obtained by changing the rotation curve of the disk. In the 
present study we have adopted a very simplified model for 
the vertical structure of the disk. More realistic models of 
the radiative transfer in the vertical direction would be 
needed to address this issue further. 

An additional difficulty with the use of the proposed 
diagnostics for mid-infrared lines is related to the fact that 
an emission line at wavelengths where the dust disk is 
optically thick (such as the lines in the infrared wave- 
length range) would be produced only if there is a spa- 
tial or temperature separation between the warm gas 
and the dust. In this respect, the recent non-detection of 
H2 rotational lines from a small sample (two Herbig Ae 
stars and one T Tau ri star) of pre-main-sequence stars by 
iRichter et al.l l)2002l) suggests that gas and dust are not 
spatially sep arated, as would resu lt in some flaring disk 
model ((Chiang fc Goldreichl IT997i) . This may indeed be 
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the case for the highly accreting FU Orionis objects that 
we consider here. 

Finally wc have come to the conclusion that optically 
thick, sub-millimetric line profiles, such as those of some 
CO lines, often observed in protostellar disks, are basically 
free from the difficulties noted for the mid-infrared lines 
and thus represent the best tool to probe the physical 
conditions of the outer parts of the massive accretion disks 
in FU Orionis objects. 
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Figure 6. Molecular hydrogen emission line profiles in the self-gravitating (solid line) case and in the strictly Keplerian 
(dotted line) case, for the conservative model with « 36 AU. The upper panel shows the S(0) line at 28/im, with a 
resolution of 5km/sec (left) and 4km/sec (right); the middle panel shows the S(l) line at 17/im, with a resolution of 
5km/sec (left) and 4km/sec (right); the lower panel shows the S(2) line at 12;tim, with a resolution of 3km/sec (left) 
and 2km/sec (right). 



14 



G. Lodato and G. Bertin: Probing the rotation curve of FU Orionis objects 




10 20 30 5 10 15 



r/AU r/AU 

Figure 7. Local emissivity as a function of radius for the S(0) (left) and S(l) (right) line profiles. Thick solid line: 
optically thin, 7 = 1; Thin hne: optically thin, 7 = 1.5; Dashed line: optically thick. The small picture in the right 
panel is a blow up of the profiles at small radii. 
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Figures. S(l) line profile in the optically thin case, with 7 = 1 (left) and 7 = 1.5 (right). 
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Figure 9. ^^CO 1-0 line profile at 110 GHz. Solid line: for a model applicable to FU Ori with a non-Keplerian 
rotation curve associated with a self-gravitating accretion disk; Dotted line: for a strictly Keplerian model of FU 
Ori. The results are shown for Vs « 18 AU (left) and Vs ~ 36 AU (right). 



